Theory of self-assembled smectic- ^4 "crenellated disks 
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Smectic- A monolayers self-assembled in aqueous solutions of chiral fd viruses and a polymer 
depletant have been shown to exhibit a variety of structures including large, flat disks and twisted 
ribbons. The virus particles twist near the edge of the structure in a direction determined by the 
chirality of the viruses. When fd viruses and their mutants of opposite chirality are mixed together 
in nearly equal amounts unusual structures referred to as "crenellated disks" can appear. These 
disks are achiral overall but the twist at the edge alternates between left- and right-handedness. To 
minimize the mismatch where the two regions of opposing twist meet, the "crenellated" structure 
exhibits cusps rising out of the plane of the monolayer. We use a phenomenological elastic theory 
previously applied to flat disks and twisted ribbons to analyze an analytic model proposed to describe 
the "crenellated" structure . When compared with flat, circular disks, we find that the model 
"crenellated disks" are stable or at least metastable in a wide region of the phase diagram spanned 
by the Gaussian curvature modulus and the edge energy modulus, with a large energy barrier 
separating the two structures. The director pattern and geometric parameters of the "crenellated 
disks" are found to be in qualitative agreement with experimental observations. 

PACS numbers: 61.30.-v, 61.30.Cz, 64.70.M- 



I. INTRODUCTION 

The smectic A (Sm-A) liquid crystalline phase is a lay- 
ered structure that expels twist and bend deformations 
analogous to the expulsion of magnetic fields in supercon- 
ductors [H,0- In a superconductor, magnetic fields pene- 
trate into the bulk phase over a distance measured by the 
London length. Similarly, in the Sm-^4 phase, two pene- 
tration depths can be defined to respectively describe the 
penetration of twist and bend deformations. In a Sm-yl 
phase composed of chiral molecules, the chirality (specif- 
ically, the chiral term in the Frank free energy) plays the 
role of the magnetic field in a superconductor. Chirality 
favors twist deformations which can appear near the edge 
of a smectic layer for low chirality. In the case of high 
chirality the twist grain boundary phase [3( appears for 
type II smectics (i.e., those with suitably large twist pen- 
etration depth). The latter phase is the smectic analog of 
the Abrikosov vortex lattice in type II superconductors. 

Single layer chiral Sm-j4 structures have drawn sig- 
nificant attention recently, motivated by experimental 
studies of membranes formed of rodlike fd virus parti- 
cles which self-assemble in the presence of a nonbinding 
polymer depletant 0, Q. Depending on the concentra- 
tion of polymer depletant the viruses self-assemble into 
a variety of structures including large flat disks (of or- 
der 10/xm in diameter), twisted ribbons and double- and 
triple-helical structures formed from ribbons. In the flat 
disks the long axes of the viruses are aligned with the 
monolayer normal in the interior of the disks. Near the 
edge of the disk the rods tilt due to both the chiral na- 
ture of the viruses as well as the tendency to minimize 
the area of the virus-polymer interface. This curved edge 
is also observed when achiral rods are studied. For achi- 
ral viruses the spontaneous twist at the edge can be ei- 
ther clockwise or anticlockwise (when the membrane is 



viewed from above), while for chiral viruses the sense of 
the twist is naturally determined by the handedness of 
the virus. Birefringence measurements yield a twist pen- 
etration depth of approximately 0.5/iin [5[. 

Recently, a new self- assembled fd structure was ex- 
perimentally discovered @ when right-handed /rf-Y21M 
viruses and left-handed wild-type fd viruses were mixed 
together in proportions such that the mixture is nearly 
achiral. Fluorescence images revealed that the two types 
of viruses mix uniformly without any phase separation. 
However, unlike the chiral case where flat membranes 
were observed, in the achiral limit the membranes are 
flat in the interior but exhibit a series of cusps along the 
edge where the surface of the membrane rises out of the 
plane of the monolayer, with adjacent cusps alternating 
above and below the plane. Between neighboring cusps 
the membrane exhibits a small radial bulge within the 
plane of the layer. These structures have been termed 
"crenellated disks" @. Images produced using the 3D- 
LC-PolScope indicate that the twist at the membrane 
edge between adjacent cusps is of uniform chirality, alter- 
nating between right- and left-handed as the edge is en- 
circled. The cusps are defects in the director field where 
right and left handed twist meet. The membrane rises 
out of the plane to reduce the director mismatch. Min- 
imizing director mismatch then requires that adjacent 
cusps alternate above and below the plane of the disk as 
shown in the schematic drawing, Fig. [1] 

Theories describing the large flat disks and twisted rib- 
bons have been constructed using the de Gennes 
model for the Sm-^4 phase generalized to include chirality, 
and in addition, in the case of the ribbons, the Helfrich 
model [ill for the surface bending energy. These the- 
ories use a simple form for the edge energy, the interac- 
tion of the rods at the edge with the polymer depletant, 
proportional to the edge length. A more realistic model 
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FIG. 1: (Color online) Schematic edge-on illustration of the 
arrangement of fd rods in the vicinity of a cusp at the edge of 
a "crenellated disk", (courtesy of Ref. 



incorporating surface tension and the "melting" of the 
smectic order at the edge has recently been developed 
19]. When applied to twisted ribbons the theory with 
the simple edge model yields good qualitative agreement 
with experimental measurements of the ribbon's pitch to 
width ratio providing that the Gaussian curvature mod- 
ulus appearing in the Helfrich energy is positive, in con- 
trast to the negative values typically measured in lipid 
monolayers or bilayers jl4| . By comparing the free en- 
ergy per unit area of a twisted ribbon with the corre- 
sponding energy of a large, fiat membrane, a first-order 
phase transition between the two structures was found in 
agreement with experimental observation. However, the 
predicted value for the edge energy modulus was found to 
be an order of magnitude lower than that measured ex- 
perimentally, presumably due to the very simple nature 
of the edge energy. 

In this paper we apply the above theory to an analytic 
model 15] proposed to describe "crenellated disks" . Be- 
cause of the complexity of the shape we do not directly 
solve the Euler-Lagrange equation obtained by minimiz- 
ing the free energy as was done in the case of twisted 
ribbons. Rather we use Monte Carlo (MC) simulations 
and determine the geometric parameters of the model 
disk and the relative stability of "crenellated" and flat, 
circular disks. We explore the stability of the "crenellated 
disk" model as a function of the Gaussian curvature and 
edge energy moduli. We find regions of this space where 
the "crenellated disk" is either stable or metastable. In 
the latter case there is a large energy barrier separat- 
ing it from the flat disk. The geometric parameters for 
the "crenellated disk" throughout the phase diagram are 
found to be in qualitative agreement with experimental 
observations, as is the director pattern of the fd viruses. 

The work is organized as follows: in the next section we 
review the theory used earlier to study twisted ribbons 
and apply it to the "crenellated disk" model. In Sec. IIIII 
we present the results of our MC analysis of the free en- 
ergy of the "crenellated disk" . We offer some concluding 
remarks in the final section. 



II. FREE ENERGY OF SM-A "CRENELLATED 
DISKS" 

Our analysis of the free energy of a Sm-^4 "crenellated 
disk" uses the achiral limit of the theory presented in 



Ref. 10] for twisted ribbons. The free energy F of 
the monolayer is the sum of the Helfrich bending energy 
[TH H3 |. the de Gennes energy Q for the Sm-A order 
and director fluctuations and a simple edge energy pro- 
portional to the edge length. Specifically, we have: 



F = J Uh + Jn)dA + 1 j 



dl 



(1) 



where fu and /„ are the Helfrich and de Gennes free 
energy densities respectively and 7 is the edge energy 
modulus (bare line tension). The Helfrich free energy 
density is given by 



f H = \k(2H) 2 + kK c 



(2) 



where H and Kg are the mean and Gaussian curvature 
of the surface respectively, k is the bending rigidity and 
k is the Gaussian curvature modulus. We have assumed 
a zero spontaneous curvature because of the up-down 
symmetry of the system. 

The achiral de Gennes free energy density in the one- 
elastic constant approximation and with the assumption 
of perfect smectic order is given by (Toj : 



/„ = ijf [(V ■ n) 2 



(V x n) 2 ] + icsin 2 6> (3) 



where 9 is the relative tilt angle of the director with re- 
spect to the local surface normal, K is the single Frank 
elastic constant and C is a tilt free energy modulus. The 
twist penetration depth is given by A t = y/KfC. 

A monolayer of general shape can be modeled mathe- 
matically as a two-dimensional surface embedded in three 
dimensions. The surface is given by a position vector 
~Y{u\,U2) parameterized by two coordinates, Ui and U2- 
To calculate the free energy of the membrane we use the 
following geometric quantities [l|| [TtJ : 

Y, = 9,Y, Y« - didjY = r^Y fc + LuN, 



diN = d N Yi = -L ij9 i k Y k , 
gij = Y i -Y j , g v =(g ij )- 1 , g = detpy, 
Uj = Ya ■ N, = {Lij)- 1 , L = detiy 



(4) 



where di = d Ui , and <9/v denotes the partial derivative in 
the normal direction N. The indices i, j, k — 1,2 and we 
sum over repeated indices. The tensors and are 
the first and second fundamental forms of the surface, 
respectively. The Christoffel symbols rk are defined by 
the relation = g km Y; L j ■ Y m . The unit normal vector 
of the surface is given by 



N 



Yi x Y 2 

V9 



The Gaussian and mean curvatures are given by: 

Kq = — , H = -g^Ln 
g T 3 



(5) 



(6) 
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and the surface area element is given by dA — ^fgdu\du 2 - formed by Yi, Y2 and N; the director n makes an angle 
The three-dimensional gradient operator is given by 9 with N. 

v ,,'-v,n, ■ m N . (7) 

After some calculation usin g Eqs. (©-(HI), the total free 
The director field n can be expressed in a local basis energy Eq. ([T]) is found to be: 



F= Y j y' >! " J ' "' VJ ' l]2 ' l^[(9ik9 j n k +g i in k T\ k )N-(2n k L jk + d j cos 6)Y,} j } v u ,1a s ,1a 2 
+ sin 2 9y / gduidu 2 + k J — y/g duidu 2 + - J (g^ L lJ ) 2 ^fg duidu 2 + 7 j 'II. 



(8) 



where djk is the antisymmetric Levi-Civita tensor. 

Henceforth, dimcnsionless units will be used where 
lengths are measured in units of the penetration depth A t 
and energies are measured in units of the Frank constant 
K. These units correspond to choosing K = C = 1 in 
Eq. ©. 

We now apply Eq. ([5]) to a model [l5[ of the "crenel- 
lated disk" structure observed in experiments. We as- 
sume that the disk lies in the x—y plane with an array of 
cusps at the monolayer edge. As noted in Sec. HI minimiz- 
ing director mismatch at the cusps requires that neigh- 
boring cusps have heights h(x,y) (in the z direction) of 
opposite signs. Likewise the alternation of the chirality 
of the twist at the edge means that cusps must appear in 
pairs. 

We model the cusps as local peaks or valleys with an 
exponential decay of the height \h\ from the center of 
the cusp [l5[. Specifically, for a single cusp centered at 
ri = (xi, yi), we assume that the height is given by 



hi(r - n) = A z cxp - 



r - r; 



(9) 



where A z is the maximum height at the peak (A z is neg- 
ative for cusps lying below the plane of the disk) , b is a 
characteristic length governing the decay of the cusp into 
the flat central portion of the disk, and r is the position 
vector in the x—y plane. For a disk with n cusps (n even) 
the height at any point r is given by the superposition 



h{v) = Y,h i {v-v i ) 



(10) 



We choose the origin of r at the center of the disk and 
the cusps are located at ri = i?(cos0i,sin0j) where R 
is the radius of the disk (excluding the bulges) and 0; is 
the polar angle location of cusp i. 

We model the in-plane radial bulges between neighbor- 
ing cusps labeled by i and i + 1 as follows 



Ri(4>) = R + A r sin 



tt(4> - 4>i) 



(11) 



FIG. 2: An example of the "crenellated disk" model with 
radius R = 5.0, A z = 1.0, b = 0.5 (see Eq. Q), and A r = 0.5 
(see Eq. {11])) in units of the penetration depth. The four 
cusps are distributed uniformly along the perimeter of the 
disk. Note that the director field is not displayed in this 
figure. 



where Ri specifies the radial coordinate of the disk edge 
between the cusps and A r is the magnitude of the pro- 
trusion. The edge of the disk is then given as a piecewise 
function of these protrusions between neighboring cusps. 
Fig. [2] shows an example of the shape of the "crenellated 
disk" model of Ref. [l3. 

With an analytic form of the shape specified, the dif- 
ferential geometry quantities defined in Eq. (£[]) can be 
computed explicitly as functions of x and y and sub- 
stituted into the free energy Eq. ([5]). In principle we 
could then follow the approach of Ref. fl0| and derive the 
corresponding Euler-Lagrange equations for the director 
field. However because of the complexity of the shape, 
the resulting equations are impossible to solve explicitly 
even using numerical solvers. Instead we discretize the 
underlying x— y plane using a square lattice of grid size 
0.05 and carry out an MC simulation at low tempera- 
ture (10 -4 in dimensionless energy units with ks = 1) 
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varying the geometrical parameters of the shape. We 
note that our dimensionless energy unit, K = 1, corre- 
sponds to approximately IOO/sb? 1 at room temperature 
using the measured value of the twist elastic constant in 
fd solutions 18j. Thus, thermal fluctuations of the direc- 
tor are negligible, as noted already in Refs. where 
the free energy was minimized and very good agreement 
was found between the predicted director pattern and 
experimental measurements. We initialized our system 
with both a random and fully aligned director field and 
found similar equilibrium states. We also carried out a 
multi-grid computation where the grid is finer at the edge 
where the directors are twisted and similar results were 
obtained. 




Experimentally [|[ the director field at the edge of the 
disk is observed to be tangent to the edge. We imposed 
this boundary condition in our model by introducing 
ghost directors along the edge of the "crenellated disk" . 
These directors are fixed tangent to the edge and interact 
with neighboring directors in the interior of the disk via 
the Frank free energy. We used a central difference al- 
gorithm to compute the derivatives of the director field. 
The total free energy was computed using numerical inte- 
gration over all the lattice sites inside the disk. Director 
defects can appear at the cusps and we include a defect 
core energy Kdefdl where 8 Z is the angle made by the 
director with respect to the z-axis. The summation runs 
over all the lattice sites whose distances from a cusp in 
both the x— and y-directions are less than one lattice 
constant. The energy K^ef is chosen to be equal to the 
Frank constant K which is unity in our system of units. 

Experimentally "crenellated disks" are observed to 
maintain constant area after their formation and thus 
as we varied the geometric parameters of the shape we 
adjusted the radius R so that the total area of the disk 
is kept constant. For each set of geometric parameters 
we carried out f 00000 MC cycles where each director is 
allowed to move once during each cycle and data are col- 
lected over the last 50000 cycles. The step size of the test 
move was chosen so that the overall acceptance ratio is 
approximately 50%. 

Since we are effectively solving the Euler-Lagrange 
equations whose solution corresponds to the zero- 
temperature ground states of the system, we collected 
the orientations of each director in each cycle and com- 
puted the averaged orientations for each of them and then 
used this averaged configuration to compute the free en- 
ergy instead of the usual method that collects the free 
energy every cycle and then averages that quantity. We 
have checked the validity of this method by considering a 
flat circular disk and obtaining good agreement between 
the results of our MC simulations and the results ob- 
tained numerically in Ref. [9| where the Euler-Lagrange 
equation was solved explicitly. 



(b) 




FIG. 3: (a)An example of the simulated director field on a 
"crenellated disk" with R = 5.41, A z = 0.6, A r = 0.4, b = 0.2 
in units of the penetration depth with four cusps distributed 
uniformly along the edge. (b)The corresponding simulated 
birefringence image where the gray scale is proportional to 
sin 2 0; 8 is the angle the director makes with the z axis. 



III. RESULTS 

We consider a "crenellated disk" with four uniformly 
distributed cusps. The area of the disk is approximately 
100 corresponding to a radius of approximately 5 in our 
dimensionless units. Significantly larger disks require 
computational resources beyond what are available to us. 
The decay length b is chosen to be 0.2, of the order of 
magnitude observed in experiments. An example of the 
simulated director field on a "crenellated" structure is 
shown in Fig. E^a). Near the center of the disk the direc- 
tors are approximately perpendicular to the x—y plane 
as expected because of the tilt energy term in the free 
energy. The boundary condition at the edge forces the 
directors there to lie tangent to the edge. It is seen from 
the figure that there is a defect in the director field at 
each cusp where the chirality of the director tilt near the 
edge changes sign from one side of the cusp to the other. 
The director twist on either side of the cusp relaxes over 
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FIG. 4: (color online) (a) Free energy per unit area of a 
"crenellated disk" as a function of height of the cusps A z 
with A r = 0.7, b — 0.2, k = 1.5, 7 = 1.6, and a uniform 
distribution of the cusps along the edge of the disk, (b) Free 
energy per unit area of the disk as a function of the size of the 
radial bulges A r with A z — 0.6, and k, b and 7 fixed as in (a), 
(c) Free energy of the disk as a function of the angular sepa- 
ration between adjacent cusps with A z = 0.6, A r — 0.7, and 
k, b and 7 fixed as in (a). The angular separation is given in 
units of 7r; e.g., a separation of 0.4 means that the four cusps 
are placed at (j> — 0.0, 0.47T, 0.87T, 1.2n. Uniform distribution 
of the cusps corresponds to A(j> = 0.5. The lines in all of the 
figures are guides to the eye. 



a length scale of order the penetration depth. Fig. [3Jb) 
shows the simulated birefringence image where the gray 
scale is proportional to sin 2 9. The bright regions corre- 
spond to the radial bulges of the membrane where the 
directors tilt more while the dark regions correspond to 
the cusps where the director is nearly vertical (i.e., point- 
ing in the z-direction) . This is qualitatively in agreement 
with experimental observations [6(. 

To understand why these achiral membranes would 
prefer to form cusps and bulges instead of remaining flat 



and circular, we analyze the contributions from different 
terms to the total free energy when A z or A r is var- 
ied. In Fig. SJa), the variation of the total free energy 
per unit area with the height of the cusps A z is shown 
with A r = 0.7 and the cusps are distributed uniformly 
along the edge. We have chosen k = 1.5 and 7 = 1.6; 
our results, however, are representative of a wide range 
of values for these parameters as we discuss below. We 
have assumed that the mean curvature modulus k is zero 
to simplify our analysis. Neither k nor k has been mea- 
sured experimentally in the fd systems but experimental 
evidence suggests that it is the Gaussian curvature that 
plays a dominant role. In the case of twisted ribbons [To| 
a small second-order Gaussian curvature term was added 
to the Helfrich energy in order to stabilize the energy. In 
the present case the negative Gaussian energy propor- 
tional to k does not lead to any instabilities and thus we 
have not included a higher order term in our analysis. 
Contributions from the different terms in the free energy 
are shown in Fig. [5Ja) . It is seen that the Frank free en- 
ergy assumes a parabolic shape like the total free energy, 
i.e., there is a minimum at a nonzero value of A z , because 
the out-of-plane cusp structure relaxes the mismatch of 
the directors in the right-handed and left-handed regions 
flanking the cusp. This optimal value of A z is not arbi- 
trarily large because a large value would lead to a large 
deformation in the director field in the neighborhood of 
the cusp. It is also seen that both the absolute value 
of Gaussian curvature (note that Gaussian curvature is 
always negative in this structure) and the length of the 
edge increase monotonically with A z . For low values of 
A z , the Frank energy dominates and at large A z , both 
the edge and Frank free energy combine to overcome the 
gain in Gaussian curvature free energy. Thus, a mini- 
mum in the total free energy is produced at A z « 0.6 for 
the the given values of k and 7. . 

The behavior of the free energy when the radial ex- 
tent, A r , of the in-plane radial bulge between neighboring 
cusps is varied is shown in Fig. 0Jb) with A z = 0.6, the 
value found to minimize the free energy. The parameters 
k and 7 were again chosen to be 1.5 and 1.6 respectively. 
Compared to the dependence of the free energy on the 
out-of-plane height, A z , of the cusps, the dependence on 
the size of the radial bulges is much weaker, possibly be- 
cause the cusps are sharply defected structures while the 
radial bulges are much smoother. The contributions to 
the free energy from the Frank, Gaussian curvature and 
edge energies are shown in Fig. [SJb) . Larger bulges in- 
troduce larger deformations in the director field far from 
the cusps because the directors are forced to follow a 
more curved edge at the boundary, but larger bulges 
can also reduce the mismatch between the directors near 
the cusps by forcing the directors on both sides of the 
cusps to become more parallel. These two competing 
effects appear to essentially cancel each other, and the 
Frank free energy is only weakly dependent on A r , with 
some numerical fluctuations. The absolute value of the 
Gaussian curvature increases with A r . The edge energy 
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FIG. 5: (color online) Contributions of the Frank (/«.), Gaussian curvature (/h) and edge (/b) energy densities to the total free 
energy per unit area for the "crenellated disk" as functions of (a) the height of the cusps A z , (b) the size of the radial bulges 
A r , and (c) the angular spacing A<f) (measured in units of n) between adjacent cusps. As in Fig. [4] we have set k = 1.5 and 
7 = 1.6. Note that the ranges of the y-axes in the three figure parts are different but the spans of energy are equal, allowing 
easy comparison of the magnitudes of the variation of each term. 



initially decreases slightly with increasing A r and then 
for A r > 0.15, the edge energy increases. The tradeoff 
between increasing edge energy and the gain of Gaus- 
sian curvature energy leads to a free energy minimum at 
A r w 0.7. 

In experiments on large "crenellated disks" cusps have 
been observed to lie on only a portion of the disk's edge 
rather than being distributed along the entire edge of 
the membrane. We consider this possibility for our small 
disks as shown in Fig. 2(c) where we plot the free en- 
ergy per unit area as a function of the angular separa- 
tion between the cusps, assuming A z = 0.6, A r = 0.7, 
b = 0.2, k = 1.5 and 7 = 1.6. The separation is plotted 
in units of it and thus a separation of 0.5 corresponds 
to the configuration where the four cusps are placed uni- 
formly along the edge. The contributions to the total 
energy per unit area from the Frank, Gaussian curvature 
and edge energies are shown in Fig. [3(c) . While squeez- 
ing the cusps into a small segment of the edge creates a 
flatter segment elsewhere which reduces the Frank free 
energy there, it also creates a more curved geometry in 
the region where the cusps are located which increases 
the Frank free energy in that location. From the general 
trend of the Frank free energy shown in the figure, we 
see that the latter effect overcomes the former so that 
the Frank free energy increases as the cusps are brought 
together. As expected, the edge length and the Gaussian 
curvature both increase in absolute value upon squeezing 
the cusps together because the membrane is more curved 
in the squeezed region and this region provides the main 



contribution to these two terms compared to the flatter 
region. For the k and 7 values selected, the trend of 
the change in the free energy is dominated by the Frank 
free energy and is minimized when the four cusps are 
placed uniformly as shown in Fig. 2(c). We are unable to 
explain the appearance of another minimum at spacing 
0.25 which may be due to the numerical fluctuations in 
the Frank free energy term. It appears that our disks are 
simply too small to exhibit cusps lying on a portion of 
the disk perimeter as seen in experiments. 

Experimental observations [6J indicate that "crenel- 
lated disks" with A z comparable to but smaller than 
the twist penetration depth should be at the very 
least metastable structures and possibly true equilibrium 
shapes under certain physical conditions. We now ad- 
dress the stability and metastability of the theoretical 
"crenellated disk" model. We consider the energetics of 
the disk as both the Gaussian curvature modulus k and 
the edge energy modulus 7 are varied. 

Our results are shown in Fig. [5] "Metastable" and 
"stable" in this figure refer to the "crenellated disk" in 
comparison with a perfectly flat, circular disk of the 
same area. Throughout the metastable and stable re- 
gions shown, A z is approximately one half of the twist 
penetration depth in accord with experiments. Our the- 
oretical prediction for A r is also approximately one half of 
the twist penetration depth while experimentally larger 
bulges have been observed. We attribute this to the fact 
that our theoretical calculations have used disks smaller 
than the real disks observed in experiments. It is rea- 
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FIG. 6: (color online)The phase diagram for the "crenellated 
disk" in the k^y space. "Metastable" and "stable" in this 
figure refer to the "crenellated disk" in comparison with a 
perfectly flat, circular disk of the same area. In the unstable 
region either A z or A r becomes very large, indicating an insta- 
bility toward other geometrical structures beyond the scope 
of our model. 



sonable to assume that value of A r will increase with the 
radius R of the disk. Thus, we argue that our theoretical 
model has produced reasonable values for both A z and 
A r compared to experiments. In the "unstable" region 
we have found that cither A z or A r becomes very large, 
indicating an instability toward other geometrical struc- 
tures which we cannot determine in our model. To get a 
sense of the energetics, consider, e.g., k = 3. For large 7 
"crenellated disks" are at best only metastable because 
the edge energy is too large to favor such a structure with 
a high edge-to-area ratio compared to flat disks. When 
7 is lowered, the "crenellated disk" can become stable. 
However, when 7 is lowered further the "crenellated disk" 
becomes unstable to other structures with even larger 
edge-to-area ratios such as helices. For fixed 7 if the 
Gaussian curvature modulus is too large structures with 
greater negative Gaussian curvature can again appear in 
place of "crenellated disks" . In most of the metastable re- 
gion, the energy barrier between the perfectly flat, round 
disks and the "crenellated disk" is about several Frank 
constants corresponding to several hundred fcgT in phys- 
ical units. 

In Ref. [lOj the transition from large, flat disks to 
twisted ribbons in the elastic theory was determined to 
occur at 7 0.28 and k w 0.1, corresponding to the very 
bottom edge of the metastable region of our phase dia- 
gram. This value of 7 is one order of magnitude lower 
than the experimentally measured value. The Gaussian 
curvature modulus has not been measured experimen- 
tally, though there is some evidence @ that suggests its 
value is larger, possibly by one order of magnitude, than 
the value found in the theory of Ref. [if]| . Thus, in our 



phase diagram we have explored values of 7 and k larger 
than those used to analyze the twisted ribbon. In the lat- 
ter analysis free boundary conditions were used for the 
orientation of the fd viruses at the edge of the monolayer 
which leads to a tilt angle 8 at the edge which is signifi- 
cantly smaller than the 90° measured experimentally. In 
the present case we imposed a 90° tilt at the edge both 
to fit the experimental observations and to incorporate a 
nonzero tilt in an achiral model which otherwise would 
have zero tilt (at least for the simple edge model used 
here and in Ref. Thus, we believe that the range 

of values for 7 and k shown in Fig. [5] are reasonable to 
explore in comparing theory with experiment. 

Given the experimental measurements to date, we have 
restricted our analysis to the possibility of "crenellated 
disks" with cusp heights A z , and radial bulges, A r , no 
larger than one penetration depth. We cannot exclude 
the possibility that "crenellated disks" with larger values 
of these parameters can be stable or metastable outside of 
the stable and metastable regions shown in the figure, but 
we have not explored this possibility due to the limitation 
of our computational resources. The "crenellated" region 
obtained by our current computations is quite possibly 
an underestimate. 

Kaplan and Meyer [n| have recently considered a 
model of an array of cusps on the edge of a flat mem- 
brane in the absence of the director field. Each cusp is 
modeled as a generic surface of revolution with negative 
Gaussian curvature (this model does not include the de- 
caying exponential of Eq. @ so as to maximize the saddle 
characteristics of the geometry) . By simultaneously opti- 
mizing the shape of the surface and its free edges Kaplan 
and Meyer obtain an analytic relation between k, k, 7 
and A ri which for k — yields: k — jA r . In the achi- 
ral limit their analysis of the energy of an array of cusps 
yields a transition to a "crenellated disk" with A r ~ 0.5 
for k ~ 1 and 7 ~ 2, in dimcnsionlcss units, of the same 
order of magnitude as our results. 



IV. CONCLUSION 

We have used an elastic theory of Sm-A monolayers to 
study a model [l5[ of the "crenellated disks" observed in 
achiral mixtures of fd viruses. The theory is the achiral 
limit of one used earlier 0, [l(| to describe flat disks 
and twisted ribbons in chiral fd monolayers. Using MC 
simulations we computed the geometric parameters of the 
"crenellated disks" , the director field of the viruses and 
the relative stability of the disks. The height of the cusps 
was found to be in good agreement with experimental ob- 
servation. Although we obtained a smaller radial protru- 
sion compared to experiments, our theoretical prediction 
is still reasonable noting the smaller size of the disks we 
have used in our calculations. The director field is also 
in good agreement with experimental birefringence mea- 
surements. The "crenellated disks" are found to be stable 
in a region where the Gaussian curvature modulus k and 
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7 are both one order of magnitude higher than the value 
previously obtained when the elastic theory was applied 
to twisted ribbons. However, the metastable region does 
include the values of the moduli used to describe ribbons. 
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